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Abstract 

A new model for stocks markets using integer values for each stock 
price is presented. In contrast with previously reported models, the 
variables used in the model are not of binary type, but of more general 
integer type. It is shown how the behavior of the noisy and fundamen- 
talists traders can be taken into account simultaneously in the time 
evolution of each stock price. The simulated time series is analysed in 
different ways order to compare parameters with those of real markets. 

Recently, several computational models trying to represent the be- 
havior of actual stocks markets have been presented. ||] ^From 
these results it seems to be a well established fact that, in order to 
obtain a good representation of the time evolution of actual markets, 
two type of traders agents must be included. On one side there are the 
so called noisy traders which are supposed to follow the local (in space 
and time) trend of the market. The noisy traders, also called followers, 
place their buy or sell orders on a given stock following the behavior 
of other (related) stocks. This kind of attitude seems to be an almost 
evident practice for someone trying to operate within a market. But, 
on the other hand, there are also fundamentalists traders which are 
considered to be responsible of the market turn offs. These kind of 
traders are supposed to know something more about the market and 
then are able to develop some kind of more sophisticated strategy to 
operate. They can take actions to buy or sell according to other indi- 
cators; the markets fundamentals. These indicators could depend on 
other type of information which usually may come from outside of the 
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system. In previous models the two types of agents are represented 
in a variety of ways. j|, ||] For instance, in reference 2 the influence 
of the fundamentalists is modeled by including a term that takes into 
account the tendency in the price changes. This effect stabilizes the 
prices. On the other hand, in reference 3 the two type of traders are 
included as separate entities. Then, it is realistic to think that both 
type of behaviors are acting at the same time and influence the way in 
which a specific stock price change. ^From the point of view of sim- 
ulations, another common characteristic of many of the models just 
presented seems to be the use of Ising like variables that represent the 
buy or sell attitude of market agents. 

Taking into account all the above mentioned characteristics, a 
somewhat different modeling approach is presented here. In principle, 
the model is based in previously reported models of opinion evolution 
in a closed community. Q But, instead of using Ising spins, here the 
community is modeled by a vector of stock prices x having N integer 
valued "Potts" components Xj, each one representing the price of a 
market asset (in arbitrary units). Then, associated with each Xi is 
a direction of movement value V{. The components V{ form a vector 
v. This components are of Ising type, i.e., they can take two values 
Vi = +1 and Vi = —1. The vectors x and v evolve in time according 
to the following dynamical rules. A randomly chosen component Xi is 
updated according to the equation 

Xi(t + At) = Xi (t) +v k (t) ; (1) 

while for the corresponding vi components the following equation is 
valid 

( v k (t) if \xi(t)\ < X th 
Vi(t + At) = l (2) 
[ - Vi (t) if \ Xi (t)\ > X th . 

The value of v k is obtained by choosing at random among the direc- 
tion values of each one of the neighbors, i.e., v k = v%—t or v k = v%+\ 
with equal probability. In principle, the algorithm described by the 
equations |l] and [2] takes into account the influence of the noisy traders 
which follow the trend of related stocks in order to buy (vi = +1) or 
to sell (vi = —1) an specific asset. However, it is not reasonable to 
think that the prices x could take arbitrary positive or negative values. 
There is no actual market following a given trend for ever. Then in 
order to take into account the influence of the fundamentalists traders, 
a threshold X t ^ is established for the absolute value of each X{. As it 
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can be seen in equation ^, if at any time \xi\ > X t h the corresponding 
direction of movement Vi is reversed, Vi — > — fj. This reversal pro- 
cedure simulates the influence of the fundamentalists traders which, 
when the absolute value of a stock reaches a value Xth, consider that 
the price is low enough so is time to buy or it is high enough and then 
is time to sell. In equations [l] and|2], At is proportional to iV -1 so the 
simulation time is incremented by one when all the stocks have, on 
average, the chance to evolve once. 

In order to analyze the behavior of the model, two representative 
indexes of the market, the mean value time series 

1 N 

XM(t) = ^2xi{t) (3) 
i=i 

and the corresponding returns or change of price, defined here as 

r{t)=x M {t)-x M {t-l) (4) 

are investigated by Monte Carlo simulations. N = 1024 and X t h = 30 
where used as typical parameters and T = 20000 Monte Carlo steps 
(MCS) where used in order to obtain most of the reported results. To 
avoid the existence of initial correlations the simulations are started 
with the components of the vector x distributed randomly between 
— 10 and 10 and with the vector v in a complete antiferromagnetic 
state. 

A typical path of the simulated price series xm{P) is shown in Fig. 
1. The sequence shows the well known noisy shape. The market 
turn offs have been indicated in the figure. This turn offs have a self- 
organized character, since they result from the time evolution of the 
model and their occurrence cannot be explicitly predicted from the 
dynamic equations. Previous analysis of actual markets time series 
suggest that the probability distribution of daily returns are not of 
Gaussian type, but they are fat-tailed. || This characteristic can be 
noticed if the probability distribution function (PDF) of the returns is 
plotted on the scale of the cumulative Gaussian distribution function. 
Normal distributions appear as straight lines in such representation 
while the fat-tail of other type of distributions result in a departure of 
the straight line. For the model presented here, the PDF of the price 
returns is plotted in Fig. 2 and the tails are clearly visible . 

According to previous results, it is known that actual markets 
time series cannot be modeled by series of independent and identically 
distributed realizations of random variables with a given distribution 



3 



(fat-tailed or not). This is a consequence of the existence of non- 
stationarity in the process that generate the series. It has been shown 
that a direct method exists in order to detect the non-stationarity, it 
is the calculation of the autocorrelation function [||, [|] 

acf (r, = ^ + M-T)(r t _ tl -T) 

of the absolute value of the returns |r(i)|. For Brownian motion the 
acf(r, t') of r(i), r 2 (t) and \r(t)\ fluctuates around zero for t' > 1. In 
Fig. 3 the acf(r, t') for the r{t) time series generated by the model is 
plotted up to t' = 50. The horizontal lines represent the 0.95 confi- 
dence interval of a Brownian random walk. Clearly, it is shown that 
the simulated series has a certain degree of non-stationarity very sim- 
ilar of those just reported for actual markets time series. 

Another parameter that can reflect the deviation of a time series 
from the Gaussian distribution is the excess kurtosis, defined as 

where /i4 is the fourth central moment and a is the standard devia- 
tion of the series under study, k is defined to be zero for a normal 
distribution, but ranges between 2 and 50 for daily returns of actual 
stock markets data have been reported. |], 0] An average value of 
TZ = 2.69 was obtained when the kurtosis is calculated on the absolute 
values of the returns. The value k was obtained by averaging over 100 
independent runs of T = 10000 MCS. 

Finally, the simulated series was analyzed using the Hurst R/S 
method. [Q] The R/S analyses begins with dividing the time series 
in segments of equal length and normalizing the data in each seg- 
ment by subtracting the sample mean. Then, the rescaled range 
(range/standard deviation) is log- log plotted against the segment size. 
By linear regression of the plot the Hurst exponent H is obtained. 
The values of the exponent reflect some characteristics of the series: 
for H > 0.5 the series is said to be persistent, if H = 0.5 the series 
represent a normal distributed random walk, while for H < 0.5 the 
series is considered to have anti-persistent characteristics. The R/S 
analysis for the price and returns series is presented in Fig. 4. The 
straight lines through the points indicates that the Hurst exponent is 
H x = 0.885 ± 0.01 for the x M {t) series and H r = 0.569 ± 0.01 for the 
r(t) series. Since the values are greater than 0.5 both series show a 
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persistent character which results from the existence of long term cor- 
relations in both distributions. In particular, the value of H r is very 
close to the value of the Hurst exponent calculated on the returns of 
the USD/DEM exchange rate. 

Although several other analysis could be made, the above pre- 
sented results show that the simulated time series obtained from the 
operation of the dynamic rules [l] and ^ can be considered as a good ap- 
proximation of the behavior of actual stock markets. Also, the discrete 
nature of the price variables follow closely the same characteristic of 
actual stocks prices. Further refinement could be made on the model 
in order to reproduce more closely some particular characteristic of a 
given market. 

This work was partially supported by a research grant from Uni- 
versidad Nacional de Mar del Plata (Mar del Plata, Argentina). 
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Figure 2: Probability distribution function of the price returns plotted on an 
inverse error function scale. The departure from the straight line show the 
existence of fat-tails in the returns distribution. 
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Figure 3: Autocorrelation function acf (r, t') (see equation f|) for the absolute 
returns \r(t)\ generated by the model plotted up to t' = 50. The horizontal 
lines represent the 0.95 confidence interval of a Brownian random walk. 
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Figure 4: R/S analysis for the price (iM) and returns (r) time series. 
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